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ABSTRACT 

I have studied condensation modes in magnetized cylindrical plasmas, 
concentrating on how magnetic field affects the stability. It is found that 
the effects of magnetic field (shear, twist, and strength) on the 

condensation modes are different depending on the wave vector. For modes 

whose wave vector is not perpendicular to magnetic field lines the plasma 
motion is mainly along the field lines; the effects of magnetic field on the 
modes are negligible except on the heat flow parallel to the field line. 
For a mode which is localized near a surface where the wave vector is 

perpendicul ar to the field line, the plasma moves perpendicular to the line 

carrying the field line into the condensed region; magnetic field affects 
the mode by building up magnetic pressure in the condensed region. 

The stability of condensation modes strongly depends on how density and 
temperature vary with field twist. The stable nature of global quiescent 
prominence magnetic configurations implies that prominences form for low 
field twist for which ideal MHD modes are stable; plasma temperature should 
increase with field twist for stable prominence formation. 

I. INTRODUCTION 

Solar prominences are a very intriguing phenomenon. They are cool 
dense material imbedded in hot and tenuous coronal plasmas. After 
Fi el d( 1965 ) various authors(Nakagawa 1970; Hildner 1971; Heyvaerts 1974; 
Chiuderi and Van Hoven 1979) have studied condensation modes in a plane slab 
geometry to understand the prominence formation. I have studied the 

condensation modes(An 1984a, 1985) assuming that prominences are formed in a 
cylindrical magnetic geometry! Anzer and Tandberg-Hanssen 1970; Pneuman 
1983). This magnetized cylindrical plasma is subject to ideal MHD as well 
as radiative( condensation) instabilities. For coronal conditions, the 
radiative time scale(tr~10 3 sec.) is much longer than the MHD time 
scale(t^l sec.), and condensation modes appear as a first order of 
£ ( = (t^,/t r ) a ~ 10“* ) (An 1985). The possible coexistance of ideal MHD 
and condensation modes with very different time scales requires a careful 
analytic manipulation for the study of condensation modes. I developed a 
mathematical techique for the study (An 1985). 

Questions about condensation modes in magnetized plasmas are how 
magnetic field (twist, shear, and strength) affects the modes and how the 
stability depends on different choice of temperature and density profiles. 
Here, magnetic shear arises from different field twist on each flux surface. 
Magnetic shear is an important stabilizing mechanism for ideal MHD 
instabilities and has been considered to be important for condensation 
modes (Chiuderi and Van Hoven 1979). I will show that magnetic structure 
(twist and shear) does not have important effects on condensation modes and 
will prove mathematically the insignificance of MHD effects on the modes for 
general magnetic field configurations. Since condensation modes are 
hydrodynamic in nature it is necessary to study the stability of an 
equilibrium with different temperature and density profiles for better 
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understanding the prominence formation and stability. 

II. DERIVATION OF CONDENSATION MODE EQUATION 

I assume that the plasma is governed by ideal MHD equation with 
radiation and heat conduction. The ambient heating rate is assumed constant 
in time and the optically thin radiative energy loss rate is used. Due to 
the two components of magnetic field (longitudinal B z and poloidal B 0 ), 
the field lines are twisted. Since MHD and radiative time scales are very 
different we make the governing equations dimensionless. Physical 

quantities, P, B, ? , T are normalized with their standard values and time t 

is normalized by the radiative time scale t r . Here radiative and MHD time 
scales are defined as t r =P 0 /R and a t 7 „=f a z /B* respectively. The 

quantity £ is defined as £ =(t w /t r ) which is much smaller than 1 for 

coronal plasmas. The quantities P r , f e , and B 0 are standard ambient 
coronal plasma pressure, density, and magnetic field and a is the radius of 
the cylinder cross section. H is the ambient heating and R is the radiative 
energy loss rate. ^ is defined as P 0 /B* . I derive the linear 

stability equation by linearizing the governing equations and assuming that 
the perturbed quantities have the form -f(r,t) = -f (r > exp {HTne* ikt) -hurt). 
The second order differential equation for linear stability has the 
following form, 


Here X is an eigenfunction and L is a second order differential operator and 
depends on w , £ , and equilibrium quantities. The condensation mode 

equation can be derived, noting £<<1, by expanding L and X in power of £ 
and by taking zero and first order equations. The explicit expression of 
zero and first order equations are as follows. 
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Prime(') denotes a derivative with r. 
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For the derivation of eq.(2) and (^) we assume that k*B 4 0 and 4*^+^^ =0 
in the plasma. The condition Tc’B40 means that the plasma should be stable 
to the ideal MHD mode. In order to study condensation modes, we have to 
solve equation (2) and (3) for a given boundary condition of X 0 and X, . 
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III. STABILITY 


The effects of magnetic fields on condensation modes vary depending on 
whether k-E is zero or nothin the plasma. Let us discuss how the magnetic 
field affects the mode of k-B=0. We cannot use Equ.(3) because it is 

derived assuming Tc*l$*0 in the plasma. According to An ( 1984b, c) the 
eigenfunction of Equ(l) becomes localized near r=r s where Tc*B=0 as m and k 
go to infinity. In this limit the growth rate of the condensation mode is 



Here, 




(v- i) 


p = PoiCl + Cf ) ; 

C 2 , = PyP 0 /Pc , Cl = Bl/p 0 ; 
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The growth rate, Equ.(4), does not have a heat conduction term because the 
term appears with "k*B, which is zero at r =r. The local mode will stay 
unstable while other global modes are stabilized by heat conduction. 
Equ.(4) shows that magnetic field directly affects the stability. If W<0, 
magnetic shear plays a stabilizing effect through Y which determines ideal 
MHD stability. As the magnetic field goes to infinity with plasma pressure 
fixed, u> becomes the isochoric mode growth rate; magnetic pressure does not 
allow plasmj condensation. 

ForK*B¥0, Equ . ( 3 ) can be used for condensation modes. I have studied 
the stability of two different equilibria with a longitudinal current 
density profile defined as ^cr) s ( l - r x )*. 

The magnitude of <* determines magnetic shear; a higher value of produces 
higher shear. Profile(A) has the density and temperature profiles defined 
as f (r>= / i^O-r^and T(r)=P(r)/ f(r) respectively while profile (B) has 
T( r ) =1— <5i(l-r A ) and f (r)=P(r)/T(r). Since plasma pressure increases with 
field twist, T(r) of (A) and P(r) of (B) increase with the twist while 
P(r) of (A) and T(r) of (B) do not change. By solving Equs (2) and (3) for 
different q we study the stability of condensation modes for profile (A) 
and (B). Here q fl is q(=2H B 2 f* /AB 0 ) at r=0 and A is an aspect ratio, q 
stands for the degree of field twist; higher q implies lower field twist. 
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Fig. 1 shows growth rate to v.s. q 0 for profile(A). Note that the 
temperature increases but density is constant as twist increases! i .e. lower 
q 0 values). Higher field twist stabilizes the condensation mode due to the 
increase of temperature (higher heat conduction and lower radiation loss). 
Growth rate increases with 5, for given q 0 because density at r=0 increases 
as 1+5, resulting in higher radiation rate. Fig. (2) shows the stability of 
profile(B). It shows that the mode becomes more unstable as twist increases 
because density of (B) increases with field twist. The growth rate of 
5* =0.53 is higher than that of £*=0.4 for a given q 0 because temperature 
at r=0 decreases with as 1- S x resulting in lower heat conduction and 
higher radiation. Fig.l and 2 have poloidal and longitudinal mode number m 
and n equal to 1, aspect ratio A=10, and longitudinal magnetic field 10 
Gauss. Because we use =0 the equilibrium does not have shear, in other 
words, every flux surface has same field twist. Fig. 3 shows the stability 
of equilibria with different values. The result seems to show that 
magnetic shear affects the condensation mode. However, a careful study 
shows that the different stability is not due to magnetic shear but due to 
different pressure profiles caused by different choice of « * . 

It is found from the resulj^that the effect of magnetic field on the 
qlobal condensation modes with k»B*0 is negligible while the local mode with 
k*B= 0 is affected by magnetic field. Why are the effects of the field on 
the two modes different? We find the answer by deriving a relation between 
parallel and perpendicular components of plasma displacements. 


«• r ® ° x + r ° ( r •* - kB^X' - Bl x B 0 ) r ~ k *' _ B ^ b /r - €Po co 2 ({ • B 0 ) . 


If k'B^O the perpendicular displacement is times the parallel 

component. Since 6 — 10~ 6 , the equation implies that plasma moves mainly 
parallel to the field line when it condenses. For 1<'T}=0 parallel motion is 
zero, in other words, plasma moves perpendicular to the field lines carrying 
the field lines into the condensed region; magnetic field affects the modes. 
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IV. CONCLUSION 


The stability of condensation modes strongly depends on how density and 
temperature vary with field twist. If plasma temperature increases with 
twist but density does not (e.g. profile(A)), then condensation modes are 
unstable for low field twist-when ideal MHD modes are stable. This result 
implies that prominences form in a globally stable magnetic configuration, 
which may explain the stable nature of prominences. On the other hand, if 
plasma density increases with twist but temperature does not (e.g. 
profile(B)), condensation modes becomes more unstable. If plasmas obey 
profile (B), we may not observe stable quiescent prominences. 

If the effect of magnetic field on condensation modes is insignificant, 
what is the role of the field on the prominence formation? We may say that 
the magnetic field has active and passive roles in the formation. The 
passive role is to insulate the prominence material from hot corona and to 
guide the plasma motion along the field lines when condensing. The active 
roles are to trap the outflowing plasmas (solar wind, spicules, and etc.) 
which then accumulate on the field lines, and to hold and support them 
against gravity. Depending on the strength and configuration, the magnetic 
field can trap the plasma effectively, enhancing the density sufficiently to 
initiate condensation, and can support the condensed plasmas to form a 
prominence(An, et.al 1986). Without the passive role, however, no 
prominence can form in a hot coronal plasma. 
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